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Abstract 

The factorial moments (FM) of the multiplicity distributions in one- 
and two-dimensions are studied for Pb+Pb collisions at 158 AGeV. Ob- 
vious intermittence slopes are observed in the relation between FM, and 
the partition number M. The Bose-Einstein correlations (BEC) are in- 
vestigated in the emulsion data sample. The influences of the two-particle 
correlation functions (TPCF) on the scaling behavior of FM are also in- 
vestigated. The dependence of TPCF on the invariant relative momentum 
Q, as well as on the pseudo-rapidity, Srj, are calculated. It is found that 
TPCF fulfill the requirements of the power-law scaling comparable with 
the scaling rule of FM. The expected TPCF contributions to FM are 
quantitatively estimated. 
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1 Introduction 

The intensity interference between the produced particles is the unique experi- 
mental tool to directly probe the space-time evolution of the interacting system. 
First, it was used in the radio-astronomy (HBT) H to estimate the astronomical 
distances. After that, the principles of HBT are used by G. Goldhaber, S. Gold- 
haber, W. Lee and A. Pais (GGLP), in order to study the angular distribution 
of produced particles in p—p annihilation |0j . GGLP describes the phenomenon, 
that the like-sign charged (identical) particles are emitted within space angles 
smaller than that of unlike-sign charged (non-identical) ones. Evidently, this 
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appears as a reason of the symmetrical properties of boson wave functions and 
therefore, the emission probability is very sensitive to the Bose-Einstein statis- 
tics. GGLP is therefore known as Bose-Einstein correlations (BEC). In this ar- 
ticle, we consider the simplest case of two-particle correlation functions (TPCF), 
namely BEC. It is known that the shape of TPCF enhancement depending on 
the relative momentum Q can be related to the emission source dimension R 
PJ and/or to it's life time r ||. In principle, the interference between identi- 
cal produced particles occurs and consequently their TPCF get powerful, when 
certain conditions are satisfactorily fulfilled. The produced particles should 
be non-coherent and simultaneously emitted, i.e. the following similar to the 
Heisenberg's uncertainty relation has to be satisfied. Ap ■ Ar < 2irh 0, ||, f|, ||. 
Ap = \pi — P2I) and Ar = \r\ — r%\ are the relative momentum and distance 
of the considered particle pair, respectively. In heavy-ion collisions the typical 
dimension of the emission source is Ar « 6 fm . Then the interference takes 
place for Ap up to 50 MeV/c [0. As mentioned above, TPCF can be used to di- 
rectly measure many important parameters (e.g. dimension, lifetime, expansion 
velocity of the emission source). These measurements are highly necessary to 
study the dynamics of strong interactions and also to probe whether the energy 
density exceeds the critical value to produce the quark-gluon plasma (QGP) 




The two-particle interferometry is recently used to study the dynamics of 
emission source until the final state of particle production [ fL0| . The Fourier 
transformation of the space-time density distribution of emission source can ex- 
perimentally be estimated by TPCF. Generally, the intensity interference is not 
a characteristic phenomenon of bosons alone. It can also be applied for the 
fermions. Of course, for fermions, the intensity correlations appear due to the 
decreased likelihood. Furthermore, the two-proton interferometry introduced 
in is an effective applicable tool at low energy. Since the protons, in con- 
trast to the pions, were present and effectively participating in all phases of 
the colliding system from the beginning until the freeze-out. The suggestion 
of the particle-antiparticle interferometry |l2| is generally based on the dis- 
crepancies between the s-channels of the particle-antiparticle, and that of the 
particle-particle interactions. 

Coming to the second topic of this article; the factorial moments (FM) which 
have found increasing interest among the high energy physicists, since the pio- 
neer works of Bailas and Peschanski during the second half of 1980's fl3|| . The 
concept "intermittence" is evidently used to describe the fluctuations in the den- 
sity distributions of particle production, i.e. the power law of FM as functions of 
the successive phase-space partitions (Jr; — > 0. To get a satisfactory explanation 
for such power-scaling behavior, many models are proposed. All these theoret- 
ical approaches are trying to give an answer to the question: "What are the 
possible sources of the power-law behavior in the multi-particle production?" 
. Following is a short list of some of these approaches. 

• Self-similar branching processes, 

• Conventional short-range correlations, 



2 



• Multi-particle cascade at the final state, 

• Critical second-order phase transition to QGP. 

Besides these dynamical sources, there are some features g for the 
responsibility of TPCF for the enhanced particle fluctuations, especially in small 
phase-space intervals. Through successive decreasing the relative momentum of 
two identical particles, an exponential increasing of their TPCF is high probably 
expected. Here we are trying to ascertain whether the reduction of the phase- 
space of pseudo-rapidity A77, which basically are performed to investigate FM, 
conducts to auxiliary enhancement of TPCF or vice versa. Experimentally, quite 
strong evidences are registered jl7j for the increment of TPCF with successive 
reduction of relative momentum Q = \p\ — P2 | , as well as with the decreasing of 
the pseudo-rapidity phase-space 6rj — \rji — 172 1. 

In this paper both subjects and their traces in Pb+Pb collisions at 158 AGeV 
are researched. For the data acquisition and the experimental setup we refer 
to [jl8| [l9| and also to the next section. Here it is important to emphasize that 
by using the nuclear emulsion one has on the one hand angular resolution and 
local accuracy as good as ~1 mrad and 1 fim, respectively. But on the other 
hand the identification of produced particles and the estimation of their charges 
and impulses can not be achieved. This article is organized as follows: Section 
U introduces the experimental setup. Section || deals with the basic definitions 
and formalism of TPCF, FM, and the strip integral correlations. Section ^ is 
devoted to the experimental results. Finally, the summary and final conclusions 
are given in Section ||. 

2 Experimental Setup 

The data used in this paper are retrieved from some of Pb chambers irradi- 
ated at CERN-SPS during 1996 for the EMU01 Collaboration. This material 
is consummately measured and analyzed at Marburg university by using our 
measuring system MIRACLE Lab (ll| [ll], |2lJ . The collisions are recorded with 
82 Pb beam accelerated to an incident momentum of 158 AGeV/c and directed 
towards stationary lead foil with thickness ~ 250 /im. The lead foil was posi- 
tioned in the front of seven plastic sheets coated on both sides by the nuclear 
emulsion FUJIET-7B. The emulsion is acting in these investigations as detec- 
tor. The emulsion has a constant area (10 x 10 cm 2 ) but different depths. The 
first two sheets have emulsions with 180 ^m in depth. The rest are coated by 
emulsions with thickness 90 /im. 

The exposure process is mainly controlled by counting the heavily ionizing 
particles by scintillator and discriminator with relatively high threshold settings. 
The number of beam particles is counted by using an additional counter installed 
behind the chamber. Scaler and driver electronics (CAEN N145) are also used 
to build up and to transfer pulses to the SPS "kicker" magnet which immediate 
consequently removes the beam when the number of particles transcends 3000 
per spot. The spot area is about 6 /urn 2 (beam density of ^5 • 10 2 nuclei/cm 2 ). 
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The number of collisions expected per spot is 17 20 . The scanning efficiency 
in the emulsion chamber is 0.75 ± 0.05 [ p0[ . 

A detailed description of the configuration and the arrangement of the emul- 
sion chamber can be taken from j2lj. Because of the relatively high transverse 
momentum of produced particles and fragmentations most of, if not all, these 
particles are expected to be concentrated within narrow forward cone. There- 
fore, their tracks can entirely be registered within the forward emulsions coating 
the other seven plastic sheets. The emulsion sensitivity for singly charged par- 
ticles with minimum ionization is found to be as good as 30 grains per 100 /im. 
Depending on the incoming momentum, the polar angle can be calculated. For 
the used data sample we get polar angle 1.3 mrad. Also, depending on the tar- 
get foil and on the topology of the microscope's field-of-view used for the image 
acquisition fll9| , the measuring system can effectively acquire secondary tracks 
with space angles 9 up to 30° (pseudo-rapidity values 77 = — lntan(6*/2) > 1.32). 
The singly charged particles are basically expected to be mixed with small con- 
tamination of e~e + pairs coming from 7-conversions in the target foil and from 
Dalitz production. As a reason of the reconstruction algorithm applied in MIR- 
ACLE Lab |l^| the tracks of these electrons cannot be grabbed. The possible 
overestimation on the observed particles density is ^2% plj. In this value all 
the contributions of 7-conversions in lead foil and in emulsion material have 
been taken into account. The efficiency of MIRACLE Lab is up to 96%. Large 
part of this 4% discrepancy is obviously coming from the tracks with relatively 
large space angle. In addition the missing measurements, frequent scattering, 
unresolved close pairs, and pair production of e~e + represent an another source 
of such discrepancy |Q. Therefore, we could consciously suggest to renounce 
the discussion of the effects of 7-conversions on FM. 

The investigation of FM is performed for Pb+Pb collisions with particle mul- 
tiplicity > 1200. From these central events only the particles observed within 
predefined pseudo-rapidity intervals are taken into account. These A 77 inter- 
vals (3 < 77 < 4, 3 < 77 < 5, 2 < 77 < 5, and 2 < ?/ < 6) are chosen at as 
well as around the mid-rapidity region of interaction. The calculation of FM is 
performed by means of Eq. @. As discussed before TPCF essentially get an 
exponential increment while decreasing the relative momenta or decreasing the 
pseudo-rapidity rather than a linear one. As discussed in |l7], ^2], TPCF are 
supposed to contribute to the intermittent behavior of FM with considerable 
values, especially for small relative momenta Q or for small phase-space parti- 
tion of pseudo-rapidity Srj. But as contemplated in Section [l the identification 
of secondary particles is highly required foremost to be able to study TPCF in 
our experimental data. The identification of produced particles by using the 
nuclear emulsion evidently is not to put into practice |l8|, [l9| |2lJ. Furthermore, 
since there was no external magnetic field applied over the emulsion chambers 
while the exposure process, measuring the particle momenta is not to be im- 
plemented too. Although all these difficulties, I have introduced in [jl4| three 
different approximated methods to adapt the calculations of TPCF to the data 
of nuclear emulsion. For this destination the three algorithms are utilized, in 
order to nearly identify the produced particles and to approximately estimate 
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their momenta. 



3 Definitions and Formalism 

3.1 Two-particle correlation functions 

By definition ||, [H], [24J, the two-particle correlation functions TPCF and their 
special case, the Bose-Einstein correlations (BEC) are given as 

C2ipuP2) = PMPM ' (1) 



C 2 (Q 2 ) = TV 



1 + Ae~ 



(2) 



A is the chaos parameter which essentially combines all effects responsible for 
the chaotic phenomena taking place in and affecting the entire dynamics of the 
emission source. For example long resonance, coherent particle production, de- 
tector resolution, contamination between non-pions and/or secondary produced 
particles, etc. Pi(pi) is the probability to detect the i-th single particle emit- 
ted with momentum pj. P2{Pi,Pj) is the probability to simultaneously detect 
the i-th and j-th two particles emitted with momentum pi and pj , respectively. 
J\f is the normalization factor depending only on the so-called "event mixing" 



(Section 3.3). By means of the covariant one- and two-particle distributions, 
the quantities Pi(pi) and P2(pi,Pj) can be expressed in dependence on the mul- 
tiplicity of one-particle N\(pi) and two-particles N2(pi,Pj), respectively pj| 

PiiPi) = £?i^- = JVi(pi), (3) 
P2(pi,p 2 ) = E 1 E 2 -^—^— = N 2 {p 1 ,p 2 ). (4) 

Then TPCF read the following expression 

n I \ N 2(pi,P2) ,_s 
C 2 {pi,p 2 ) = , -. , r - (5) 
Ni{pi)Ni(p2) 

The experimental steps to estimate TPCF and afterwards to calculate the pa- 
rameters J\f, A, R and/or r can briefly be summarized as follows: 

1 . Empirical estimation of C 2 by using the particle multiplicity. This esti- 
mation is simply the ratio of the correlated to the non-correlated particles 
as given in Eq. (||). 

2. Choosing an analytical expression for C2, for instance Gaussian distri- 
bution which can be given in dependence on the dimension of emission 
source R (Eq. (|J)) or also on the life time r. Depending on the modeling 
of the emission source density distribution we can extend these analytical 
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expressions for different dimensions to retrieve estimation for other pa- 
rameters. For the emulsion data I have driven a special expression given 
in Eq. © Q. 

3. Now the last step is more or less clear. Fitting the experimental results of 
Ci (step #1) with the analytical ones (step #2), to directly get estimation 
for Af, R, t, etc. 



3.2 Factorial moments and intermittence exponents 

If the phase-space for instance the pseudo-rapidity Ary is split into M bins of 
equal sizes Sr/ — Arj/M and if n& being the multiplicity in fc-th bin, then the 
"vertical" FM are given as Q 

<F q (M) >= <^K-l)---K- g +l)> 

< rik > is the average multiplicity of the A:-th bins of all events. The self-similar 
density fluctuations for the successive partitions leads to the following power-law 
scaling. 

<F„(M)><xM+<. (7) 

The parameters <\> q which represent the slopes of the relations logF g vs. logM 
are called "intermittence exponents" . Their relations with the anomalous di- 
mensions ]2q| are only depending on the orders q 

<, - jij. ( 8) 

Therefore, for multi-fractal processes d q can simply be given as functions of 
q alone. d q = di characterizes the mono-fractal processes, i.e. the particle 
production can be characterized by one fractal dimension. 

3.3 Strip integral correlations 

The integrals of g-particle inclusive rapidity densities p q {r]i, 772, • • • ,f]q) are con- 
ditional to FM given in Eq. @ Q || 

/ 5,1 \ 
p / I dmdV2 ■■■drjg p q (771, 772, ■■■,y q ) \ (Q] 

q \ <n>i J ' [ ' 

For the ideal case that there is no correlation between produced particles p q (771 , 772 , ■ • • , r) q ) 
pi(?/i)pi(??2) • • ■ Pi(Vq) and consequently, < F q >^ 1 

<FJ8r))> = L (10) 

qy " NORM v ' 
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NORM =< J Q YU drji p\{r}i) Pi(f7a) ■ ■ ■ PiiVq) > is estimated by the so-called 
"event mixing" [^7|. Generally, the estimation of "NORM" plays a very im- 
portant role with respect to the final estimation of TPCF. In-offcnsively, one 
would suggest to replace "NORM" in last equation by the distribution of non- 
correlated particles. The difficulties to determine both correlated and non- 
correlated particles in the emulsion data are frequently discussed before. "NORM" 
is our data sample — non-correlated particle-pairs — is estimated by using the 



three methods given in Q (Section 4.2.2). < J Q ■ ■ • > in Eq. (M) represents 



the correlated particle-pairs. 

Using the strip integral correlations (Eq. (|lo|)) instead of the conventional 
sum over the bin multiplicity (Eq. (^)) enables us to entirely avoid the large 
disadvantage known for the second one, namely the artificial fake partitions 
[p3| . In addition to these, the strip integral correlations have other important 
advantages to decrease the statistical errors and to increase the correlation size, 
etc. |27j . To be able to calculate the integral correlation functions in the pseudo- 
rapidity phase-space At] obviously, we would need to extend the strip integration 
to another domain, individually depending on the pseudo-rapidity variable r\ 

( / II drhPgivi, 1 ^, ■■■Vq)) 

< C ^"> >= — NORM " <"» 

Besides the extension into 77-domain performance the counting of g-tuplesQ 
in last equation represents the other difficulty. To avoid this, we use the 
Grassberger-Hentschel-Procaccia integrals |27], ^9|. Therefore, the last strip 
integration for instance for the second order (q — 2) reads the following form 



(2! e n « («»-k. - %V I) ) 
<^"»- V '""""norm L W 

n is the total particle multiplicity inside the whole considered interval Ary and 



Ti. is the "Heaviside" unit step function. Analogue to Eq. (12) and by utilizing 



suitable algorithm to count the strip integral correlations it is possible to deduce 

lr The counting algorithm used in this article has the following instructional iterations 

• Determining the phase-space Ar] and the partition width Stj. Within each phase-space 
partition (Srj), the differences between the pseudo-rapidity values r]ij = Jji — i < j 
are calculated and then compared with the pre-specified width 8rj according to the 
inner product if Eq. ( |l3[ ) . 

• Counting the q-tuple of {riij 1 , r]ij 2 , • ■ • , r\ij q } from each i-th event of the used data 
sample by using the Heaviside unit step function. 

• The results of last counting is then multiplied by q\ to take into account the number of 
permutations within each q-tuple. 

• Repeating the last steps for all available (™ ) — / tuples and then calculating the average 
over the total events. 
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a general expression for any order q > 2. This expression can also be given in 
dependence on the relative pseudo-repidity 8-q (Eq. (fjjl)) or on the relative 
momentum Q (Eq. (|l5|)). For the direct relation between both variables Q and 
Srj we refer to Note #c 



of this relation, we are apparently succeeded to derive an expression to study 
the contributions of TPCF to the power-law scaling of the conventional FM. 

3.4 TPCF as a possible source of intermittence 

Once more the conditions of the validity of two-particle interferometry are given 
in Section |l|. Since the breadth of the emission source is limited within the finite 
source diameter 2R, then TPCF (BEC) are effective, if the relative momentum 
fulfill the following relation 



In addition, from last expression the particle interferometry (non-coherence phe- 
nomenon) strongly takes place, if the breadth of the emission sources R — > 0. 
Of course we still remember the properties of FM which are suggested as an ef- 
fective signature for the non-statistical fluctuations at the final state of particle 
production |l3| and then as a diagnostic tool to confirm the QGP formation. 
FM show power-like behavior for the consecutive dividing of the phase-space 
Ar] into 5r] — > sub-intervals. Equation (|l4|) describes almost the same behav- 
ior, namely while the momentum difference is successively decreased Ap — » 0, 
an exponential increasing of TPCF is strongly expected. This behavior remains 
valid, even if the breadth of emission source R has a constant value (fixed source 
diameter). In case that the dimension of emission source sanctimoniously suffers 
a successive explosion or shrinking R — > it is obviously expected that TPCF 
become more stronger. Practically, the power-like behavior of two-particle cor- 
relation functions C2 has been investigated |l7j with the successive reduction of 
the momentum difference Ap = Q — > as well as with the successive reduction 
of the pseudo-rapidity interval Sr/ — » 0. This leads to the following power law 
which apparently is to be compared with the scaling rule given in Eq. (Q). 



TPCF are, as their name announces, the correlations between the particle-pairs 
in a very restricted region (Ar < 6 fm and Ap < 50 MeV/c). Therefore, it is not 
more difficult to explain why they are supposed to contribute to the compulsory 
correlations between the produced particles? Especially, if the investigation is 





(14) 



C q (AQ) oc (AQ) 



(15) 
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being performed within the region defined above, in which TPCF are effective. 
For all these reasons it is now more or less clear, if we predict that TPCF can 
be suggested as possible sources of the "intermittent behavior" in heavy-ion 
collisions |j§ 

4 Experimental Results 

4.1 Results of factorial moments 

In all evens of the data sample the four pseudo-rapidity intervals Ar/ defined 
in Section ^| are successively split into M — > 00 sub-intervals (bins) with equal 
widths Srj — Ar]/M — ► 0. In each such bins the particle multiplicity is counted 
and then the g-order FM are calculated according to Eq. (Q) . Figure [l] shows 
\ogF q as functions of logM^ for the orders q = {2,3,4}. The subscript in M v 
means that the phase-space partition has been performed in Ty-dimension only. 
Obviously, the experimental results can good be fitted as straight lines according 
to the scaling rule Eq. (^). The slopes represent the so-called intermittence 
slopes 4> q . The four pictures show obvious positive linear tendency. Also, it 
is to realize that the values of F q (M rj ) and </> q increase with decreasing the A 77 
intervals. The largest values are obtained within the smallest interval 3 < i] < 4. 
There is another remarkable finding that (j> q increase with increasing the orders 
q (review Eq. d)). 

By continuous increasing the partition number M v , FM exponentially grow 
(Fig. ||). The reasons for this fast upwards increasing are discussed in details 
elsewhere |23|]. Here the discussion is restricted only on the phenomenologi- 
cal behavior of the increasing of F q (M v ) with increasing the partition number 
M„ . Naturally, increasing Ar/ intervals enables us to divide the one-dimensional 
phase-space (77) into large partition number M^. Once again we also notice in 
this figure that the range of the linear dependence becomes larger with wider Arj 
windows (pictures from top to bottom). This linear dependence are studied in 
the previous figure. In the top most picture, the curvature upwards increasing 
has sharp spikes. Plateau begins to appear with increasing the Arj intervals. 
In the largest interval (bottom most picture) we have the largest plateau cov- 
ering the distance from M v « 200 to M n w 600. Beyond this saturated regions 
F q (M r) ) continue their exponential increasing. This time steeper than before. 
We continued the partition process S-q — > as long as FM increase. 

4.2 Results of TPCF 

4.2.1 Adapting TPCF in emulsion data 

As shortly introduced before it is expected that TPCF contribute to FM with 
considerable values, especially, if FM are studied within small phase-space inter- 
vals where TPCF are expected to be strongly enhanced (Section [j]). Therefore, 
before the announce whether non-statistical fluctuations (e.g. dynamical) are 
observed as the case in Fig. and §, one has to study all possible sources which 
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Figure 1: The factorial moments F q for the orders q = {2, 3,4} 
are plotted as functions of small partition number M v . The data 
used here are Pb+Pb collisions at 158 AGeV with total multi- 
plicity > 1200. The results can good be fitted as straight lines 
(Eq. ffp). Their slopes represent the "intermittence exponents" 
(p q . Although the ordinates are also log scales, their labels are 
given. 

may be responsible for such abnormal scaling behavior. After that if the fur- 
ther scaling of FM still establish linear positive tendencies with the partition 
number, this might be then referred to self-similar branching processes and/or 
to critical second-order phase transition to QGP |2^, j3l}] . Regard being had 
to the nature of TPCF alone, one can expect an exponential dependence with 
the decreasing of the considered phase-space rather than a linear one (review 
Eq. (^)). As mentioned in Section |l] the identification of produced particles 
and the estimation of their momenta are highly necessary to be able to study 
the two-particle intcrfcrometry. The two-particle intcrferometry known as BEC 
(Section represents the simplest case of TPCF. The identification of produced 
particle by using nuclear emulsion obviously is not to be put into practice. In 
three different methods have been introduced and examined to make it pos- 
sible to estimate TPCF for the emulsion data sample. For these algorithms the 
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Figure 2: TVie factorial moments F q for the orders q = {2, 3, 4} 
are depicted in log-log charts as functions of large partition num- 
ber M v . The partition process is repeated as long as FM values 
increase. An exponential growth of FM is noticed for relatively 
large M v . Within the largest An interval, we notice that FM are 
saturated. 



expression given in Eqs (|[) and @ must be modified as follows: 



1 + N 1 + Ae — ~ 



(16) 



The prediction of this approximation are compared with the ratios between cor- 
related and non-correlated multiplicity of particle-pairs (Eq. (||)). The super- 
script in Eq. (16) says that the investigation and consequently, this analytical 
form are applied for emulsion data sample. The subscript refers to the two- 
particle correlations BEC. This will be replaced by q, when we will modify Eq. 
( p^[ ) to be compared with FM for the different orders q. The results of Eq. ( fl6| ) 
will be depicted in Fig. |[ The fitting parameters TV, R, and A are compared 
with the parameters of data samples acquired by other detectors which are in 
the position of being able to identify the produced particles and to measure 
their impulses p2|. 
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4.2.2 TPCF as function of Q 

Figure || shows the dependence of TPCF on the invariant relative momentum 
Q of all possible particle-pairs of our Pb+Pb data sample observed within 
the pseudo-rapidity interval 1 < r\ < 10. For this data set we counted about 
~3 • 10 6 particle-pairs. The averaged transverse momentum < p t > is assumed 
to be constant ~350 MeV/c. By using this constant value and the related 
equations given in |Q it is possible to assign to each particle a stochastically 
calculated momentum. Then it is possible to calculate the relative momentum 
Q for each particle pair. The algorithms discussed in are used to count the 
non-correlated particle-pairs and then to classify their multiplicities according 
to Q. Counting the correlated particle pairs is now more or less unsophisticated 
process. In fact the algorithm which randomizes the angular characters of the 
original events, is preferentially used. Obviously, the randomized particles in 
this regard represent the background distribution "NORM" . In such a way the 
normalization factor TV gets an average value fixed around the unity. It is 
necessary to point out here that no correction for the Coulomb repulsion has 
been considered, since the amount of Gamov correction factor |53| is ^10% 
for Q w 1 MeV/c. With increasing the values of Q Gamov correction factor 
exponentially decays. As discussed in Section |l| the two-particle interferometry 
must be taken into account only for Q < 50 MeV/c. Therefore, the ratios of 
correlated particle pairs (multiplicity of original events) to the non-correlated 
ones (multiplicity of background distributions) for Q > 70 MeV/c can be fixed 
to the unity, i.e. for Q > 70 MeV/c, we can assume that there is no difference 
between the multiplicity of correlated and non-correlated particle-pairs. 
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Figure 3: C|' m ((5) are depicted as functions of the invariant 
relative momentum Q for the interval 1 < r\ < 10. The fitting 
parameters obtained by using the analytical relation Eq. are 
given at the upper right corner. 
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The experimental results shown in Fig. |^ are fitted by using Eq. (16). The 
obtained fitting parameters are R = 6.686 ± 0.607 fm, A = 0.428 ± 0.032, and 
N = 1.010 ± 0.007. All these values are comparable with the results obtained 
from data sample of detectors able to identify the produced particles and to es- 
timate their charges and momenta ]32j . For our one-dimensional analysis the re- 
sulting parameters are assumed to combine all their corresponding components. 
Using one-dimensional analytical expression for the correlation functions C as 
in Eq. ( |l6| ) which depends only on the invariant relative momentum Q obviously 
does not enable us to directly investigate the behavior of the other components 
of R and A with Q. The breadth of emission source R = 6.686 ± 0.607 fm is 
on the one hand evidently comparable with the r.m.s. nuclear radius of 82 Pb, 
7.11 fm on the other hand it shows that the method used here to calculate BEC 
in emulsion data sample has a very good physical acceptance. The value of 
the chaos parameter A < 1 is an evidence for the partial non-coherent emiss ion 
source and for the simultaneously emitted particles. As given in Section |3.l| 
this parameter includes all possible chaotic effects taken place in the emission 
source. More details about the results of BEC in the emulsion data sample can 
be obtained from [[hi] . 

4.2.3 TPCF as function of r) 

The conventional TPCF dependence on the invariant relative momentum Q 
does not enable us to directly^ study their contributions to the intermittent 
behavior observed in the experimental data (Fig. |l| and ||) . For this purpose, it 
is needed to study as evident as possible the direct TPCF dependence on the 
pseudo-rapidity r\ or alternatively to investigate their behavior within the region, 
Q < 50 MeV/c. Such a way it will be possible to check whether their behavior 
can be given by a power law similar to that of FM, i.e. whether TPCF are also 



intermittent. Section 4.2.4 is devoted to the second suggestion. The behavior 



for large interval 0.5 < Q < 100 McV/c will be studied and compared with the 



results given in Section 4.1. In this section we study the direct dependence of 
Of m on the pseudo-rapidity phase-space Srj. 

Figure ||A shows the dependence of Of™ 1 on 5r) for the pseudo-rapidity in- 
terval 1 < rj < 10. The intervals between the successive iterations of <5r;'s are 
constant (0.1). Taking the different statistical acceptances into account we no- 
tice that the step- wise successive increasing of Si] = — rj 2 \ leads to slow linear 
decreasing of Gf m until the value Sr] w 2.0. This value is corresponding to space 
angles 9 w 15°, i.e. TPCF of particle-pairs distributed up to ^15° are able to 

2 Since the invariant relative momentum Q on the one hand depends on the squared 
invariant mass Q 2 = Mf - (im) 2 , it can be used to study the fluctuations due to cluster 
resonance decay. On the other hand this relation reads 

Q 2 = bit + Pit + l-PltPlt cos(<5</>)] 

- [M 2 t + M| t + 2 M 2 t Ml cos (8 V )] . 

This relation implies that for given 8<f> = |0i — <fe|, Q = \pi — P2I can be given in dependence 
on Sr] = \rji — 772 1 - Such a way TPCF can also be studied in dependence on the phase-space 
partition of pseudo-rapidity 8r]. 
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V 2 < Tl < 5 
□ 3 < Tl < 5 
• 3<T|<4 



0,1 1 
Mr TUl 

Figure 4: 27ie dependencies ofC^dvi — Vs]) on the n differ- 
ences \rji — 772 1 are drawn in A. C^ m have values greater than 
1 for the r\ differences < 2.0. After that their values decrease 
exponentially. The same relations are drawn in B for the other 
Ar/ intervals. The linear dependencies are also still to be noticed 
here. Also the exponential decrease for Srj > 2.0. 



show a power-like scaling. Almost the same behavior we observed in Fig. [lj. It 
is clear that £77= \t]i — 772 1 used here does not directly coincident with the bin 
size Si] = Ar]/M used in Fig. ^and ^|. This linear decreasing with increasing 8r] 
obviously characterizes the expected scaling rule. Therefore, it can be compared 
with that of the FM (Eq. ([?])). As given before the intermittence slopes <pq can 
be directly retrieved from such log-log scale. If we continue to increase Srj, the 
decrease of C® m becomes exponential. And TPCF become smaller than unity.P 
For very large Sij some kind of plateau (saturated regions) are formed. In Fig. | 
we showed that TPCF or BEC are applicable only for Q up to ~50 MeV/c. 

3 The ratios between the correlated and the non-correlated particle— pairs define the two- 
particle correlation functions C2- For the value C2 < 1 we can suppose that the ratios and 
the corresponding functions are not valid any more. 
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Figure |I]A can be used to visualize the conclusion that the functions C^iSn) 
can be compared with C|' m (AQ). We can therefore conclude that the func- 
tions C(5rj) are to be used for Sn < 2.0 and eventually, for AQ < 50 MeV/c. 
That the functions C2{Srf) have linear dependence on the Sr) differences until 
Sr) ~ 2.0 leads to the assumption that until this value we can describe their 
dependencies by using a power-law scaling. Then the two values Sr) rj 2.0 and 
AQ ks 50 MeV/c are coincident with each other. They apparently represent the 
upper limit to consider the TPCF effects. Beyond these values TPCF do not 
affect the observed FM. 

The bottom picture (Fig. ||B) depicts the same relations with the same step- 
width 0.1 but for the other An intervals. The linear dependencies depending on 
the Arj intervals^] are also still to be noticed here. They expand also until the 
value Sr] r* 2.0. Although, C2 m {Sr]) for the small Arj intervals leave the linear 
dependence earlier than for the large ones evidently, almost the same behavior 
is also noticed here. There is an interesting finding that the increasing of Sr/ 
differences results a giant decreasing of the correlation function C|' m (<5?7). In 
this regard we can conclude that the effects of TPCF exponentially weakened 
with increasing the differences Sr]. Especially, within the two large intervals 
2 < n < 5 and 2 < r\ < 6 we have similar saturated regions as in the top picture. 
So far we can summarize that the comparison between this figure and Fig. [j] and 
U shows that our investigation of TPCF in dependence on Sr] differences results 
essential information about their contributions to the intermittent behavior. 



4.2.4 Strip integral correlations 

Using the strip integral correlations given in Section || we can calculate FM 
in dependence on the invariant relative momentum Q. This represents an ex- 
tension of the second alternative method suggested in Section 4.2.3 to check 
whether the behavior of TPCF for the orders q = {2,3,4,5,6,7} can be given 
by scaling rules similar to that of FM. The results of strip integral correlations 
are given in Fig. 4.2.4A for the interval 2 < rj < 6. Starting with a certain 
momentum difference we compare it with the calculated relative momentum 
\Pi - Vi I by means of the Heaviside unit step function H as given in Eq. @. 
We can clearly realize that the increasing of the momentum differences, i.e. 
increasing the distances between the particles, within which we are searching 
for the possible particle-pairs which fulfill the condition Q 2 — {\pi — P2\ 2 ) > 
results an additional decreasing of FM. Taking into account the effects of the 
inner product in Eq. (|l3|) the proportional increasing of FM according to the 
orders q can be understood. 



In light of the scaling Eq. ( |15|) we study the dependence of C^ m {AQ) on AQ 
for relative momenta ranging from 0.5 until 100 MeV/c. Obviously, we notice 
that the ability of Eq. (|l5|) to describe these results is limited. Therefore, it 
should be modified to simulate these non-linear relations. The validity of Eq. 



4 This generally valid for all Arj intervals. For the first two large intervals as well as for 
the interval 1 < rj < 10, evidently this value is included. For the other two intervals it is 
excluded, since the differences \rji — T72 1 are not enough to reach this value. 
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Figure 5: For the orders q = {2,3,4,5,6,7} the strip den- 
sity integrals C^ m are drawn in A as functions of the distances 
\pi — P2 1 ■ The experimental results are fitted by using scaling 
rule given by Eq. The fitting parameters are given in 

Table ^|. The same relations are depicted in B for the region 
1 < \p\ — pal < 50 MeV/c. These can good be fitted as straight 
lines (Eq. fttdj)). Their fitting parameters are listed in Table 
The relations between the q-order C^ m and the second-order 
ones are illustrated in C. The results are fitted according to Eq. 

The fitting parameters are listed in Table |^. Although we 
also use here a log-log scale, the labels of the abscissa are given. 



L5|) will be discussed later. We would therefore suggest the following empirical 
power law 



C* m (AQ) 



b ■ (AQY 



(17) 



Using this scaling rule the experimental results in Fig. 4.2.4 A are fitted and the 
resulted parameters given in Table 0. We generally notice that for all relations 
the fitting parameters increase with increasing the orders q. There is only one 
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exception for q — 2, c is too large meanwhile b is too small. By means of this 
scaling rule it is not possible to retrieve any information about the contribution 
of TPCF to the observed FM, since the intermittence slopes are entirely hidden. 
For this destination we should first divide the region 0.5 < AQ < 100 MeV/c 
into sub-intervals possibly with common properties. As discussed before TPCF 
have a restricted region outside of it they are not applicable. Last figure serves to 
view these properties. Visually, one can notice for instance for AQ > 50 MeV/c 
the values of log C® m are constant for all q. 

Although we could not retrieve more experimental points for the region 
AQ < 1 MeV/c, it is to realize that within this region C^ m (AQ) exponentially 
decrease with increasing the AQ differences. Although this decrease continues 
even until AQ ps 50 MeV/c, we can notice that the decrease within the region 
1 < AQ < 50 MeV/c is not exponential. We still remember the properties of 
the functions C^™ within this rest ricte d region. They are comparable with the 



factorial moments given in Section 3.2. As we noticed they showed an obvious 



intermittent behavior. Therefore, the subscript q gives the order of FM. For 
AQ > 50 MeV/c we can assure that Cf m (AQ) w 1.0 for all AQ, i.e. TPCF 
have no more effects on the observed particle correlations, if the two-particle 
relative momenta are greater than 50 MeV/c. These facts are visualized in this 
figure. We obviously notice that all experimental point are saturated. 

To study the behavior within the region AQ < 1 the experimental data 
should be first percolated from the Coulomb repulsion. The Gamov correction 
factor which is usually used to measure the possible interactions between non- 
coherently emitted particles, i.e. particles simultaneously emitted with almost 
the same momenta, is about 10% for AQ w 1 MeV/c. It decreases exponen- 
tially with increasing the relative momenta AQ. It has a negligible value for 
few MeV/c. The experimental results used here are not liberated from these 
interactions. Therefore, we can exclude this small region while studying the 
dependence of Cg Tn on Q. 



Table 1: The relations between \ogC(AQ) and log AQ drawn in Fig. 4-2. 4 A for 0.5 < 
AQ < 100 MeV/c are fitted according to Eg. ftP\). The fitting parameters a, b, and c 
are listed here. We notice that all values increase with increasing the orders q. 

q a b c 

2 1.188 ±0.002 0.022 ± 0.006 2.882 ±0.984 

3 1.225 ±0.066 1.654 ± 0.145 1.004 ±0.116 

4 1.365 ±0.136 3.263 ± 0.339 1.335 ±0.144 

5 1.496 ±0.211 4.938 ± 0.546 1.607 ±0.156 

6 1.606 ±0.277 6.405 ± 0.723 1.721 ±0.161 

7 1.721 ±0.349 7.412 ± 0.913 1.792 ±0.175 



Thence, we are left with the region 1 < AQ < 50 MeV/c, only. For this re- 
gion the dependence of Cg m (\pi — P2\) on the relative momenta AQ = \p% — _p2 
are drawn in Fig. 4.2.4 B. We notice that the relations in this log-log scale can 
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Table 2: Wi thin the region 1 < AQ < 50 the relations \ogC(AQ) vs. log AQ depicted 
in Fig. (.2.4B can good be fitted by Eq. m_a). We notice that the intermittence slopes 
(f> q increase with the increasing orders q. The constants k represent the intersects. 





k 


0q 


2 


0,083 ±0,001 


0,009 ±0,001 


3 


0,426 ±0,018 


0,241 ±0,020 


4 


0,631 ±0,024 


0,372 ±0,026 


■5 


0,776 ±0,026 


0,461 ±0,026 


6 


0,874 ±0,026 


0,514 ±0,028 


7 


0,937 ±0,020 


0,539 ±0,022 



good be fitted as straight lines (Eq. (jlq)). This behavior is entirely comparable 
with all pictures in Fig. |l| and in the linear regions in Fig. |]. As given before 
the tendencies represent the intermittence slopes (j) q which are listed in Table 
The increasing the orders q the increasing the slopes <f> q . This behavior char- 
acterizes the intermittent interacting system. Within this region we can then 
quantitatively estimate the contributions of TPCF (BEC) to the corresponding 
FM. To do this we should first study the dependence of the ratios 4> q /4>2 on the 
orders q. Motivated by the Ochs-Wosiek log-log scaling p3| by the adoption of 



C™_to the conventional FM given in Section |3.2| and by the limited validity o:' 



Eq. ( |15| ) we would therefore recommend the following power-law scaling to be 
applied within the region 1 < AQ < 50 MeV/c 

logCf"(AQ) = h. . log cf "(AQ) + k . (18) 



Table 3: The dependencies of q-order m (AQ) on the second-order ones de- 
fer the restricted region 1 < AQ < 50 MeV/c arc fitted 
ij). We clearly notice that the ratios r q = 4> q /4>2 increase with 



picted in Fig. 4.2^ 
according to Eq. (|L 

increasing the orders q. Defining the variables m q — r q — (|) and n q = r q — (q— 1) 
we compared them with r q . Also the differences between the successive orders 
of m and n are given. 



q 




-k 


m q 


m q - m q -i 


n q 




3 


27, 52 ±2, 18 


1,85 ±0,17 


24 


23 


25 


24 


4 


42, 37 ±2, 99 


2, 88 ±0,22 


36 


12 


39 


14 


5 


52, 43 ±3, 63 


3, 56 ±0,27 


42 


6 


48 


9 


6 


58, 43 ±4, 03 


3, 96 ±0,30 


43 


1 


53 


5 


7 


61, 17 ±4, 33 


4, 13 ±0,33 


40 




55 


2 
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The relations between the g-order Cf m {AQ) and the second-order ones are 



given in Fig. 4.2.4C. As expected all these log-log relations have positive linear 
dependencies where the ratios of intermittence slopes 4> q /<p2 can directly be 
retrieved from these relations. 4> q /4>2 provide important information about the 
reaction dynamics and also about the properties of intermittent behavior pl[ 
p3| , pl| . The fitting parameters according to the scaling rule Eq. (|l^) are listed 
in Table ||. We notice that 4> q /<f>2 apparently increase with increasing the orders 
q. These ratios are too large in comparison with (|) which are related to the self- 
similar process. They are also larger than the values resulting from the scaling 
rule of critical phenomenon characterized by (q — 1). For the same data sample 
we could proof on the one hand that the intermittence slopes <f> q / '4>2 cannot be 



simulated as self-similar process [|23, 31 . Their values are smaller than (2). On 



the other hand we could partially fit the experimental results according to the 
critical phenomenon (q — 1). In light of these results wc try in next section to 
give an answer to the question "How much can TPCF contribute to the observed 
FM?". We will measure their contributions to the most important parameter 
of the intermittent behavior, namely, the ratios 4> q /4>2- 

4.2.5 TPCF contributions to FM 

We define r q = 4> q /4>2, m q — r q — (|), and n q — r q — (q — 1). The differences 
between the successive orders of m and n are also calculated and listed in Ta- 
ble [| Taking into account the assumption that the contributions of TPCF are 
linearly added to the observed FM. Furthermore, we know that the experimen- 
tal values of 4> q / 4>2\exp obtained from FM are calculated over wide ranges of 



pseudo-rapidity phase-space Arj (Section hi ) . Then we can suggest to extrapo- 
late them into the relatively small restricted region of the valid TPCF (Srj < 2.0 
or AQ < 50 MeV/c). Therefore, we can suppose that the values of 4> q /4>2 |exp are 
still valid within the region 1 < AQ < 50 MeV/c. In the language of pseudo- 
rapidity we suppose the validity of </> g /</>2|exp for Sr] < 2.0. As given above both 
(I), and (q — 1) could be used on the one hand to simulated the experimental 
results of FM (<fi q / 4>2\exp) ■ On the other hand both of them are unable to fit 
the results of r q . We noticed from Table ||that r q are larger than both (f) and 
[q — 1). According to these assumptions we can suppose to subtract 4> q /4>2\e^p 
from r q observed within the region of effective TPCF. Such a way we get a 
primitive estimation of TPCF contributions to FM. The increasing the orders 
q the increasing both m q and n q . For (|) and (q — 1) we can realize that the 
differences between the successive orders increase with increasing the orders q. 
In contrast we notice for r q and correspondingly for both m q and n q that the 
differences of their successive orders decrease with increasing q. For large orders 
we expect that m q — m q ^\ and n q = n q -\. This can be referred to the limited 
rapidity values of the used data sample and to the restricted available dimension 
of the emission source. 
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5 Summary and Final Conclusions 



In this article we have studied the non-statistical fluctuations in Pb+Pb colli- 
sions at 158 AGeV. First, we would like to sum up the results obtained so far. 
The analysis of FM in one-dimension shows that the interacting system is ob- 
viously intermittent. The phase-space partitions are performed as long as FM 
increase. For large partition number the increase of FM is exponentially binding 
upwards. Part of this fast increasing can be understood as TPCF effects. For 
large A 77 intervals (wide data sample), saturated regions are formed. By contin- 
uous increasing the partition number FM continue to rise fast. The difficulties 
to apply BEC in the emulsion data are shortly discussed. Nevertheless we ap- 
plied one of the methods suggested in |H . Using a special analytical expression 
for TPCF we could retrieve estimations for the dimension of emission source 
R and for the chaos parameter A. As expected the randomization parameter 
A w 1. This values is an evidence for partially non-coherent emission source 
and for the simultaneously emitted particles. For the dimension R the value we 
get is comparable with the r.m.s. 82 Pb radius. These results further the use of 
our method to estimate BEC in the emulsion data sample. This enabled us to 
directly pass to the main destination of this article. The study of the TPCF con- 
tribution to the observed FM. First we studied the dependence of TPCF on the 
phase-space of pseudo-rapidity (on the differences, 5rj). For small Srj differences 
we get linear dependencies for all A77 intervals until Sr] ss 2.0. After that the 
dependence becomes exponential. This limitation is coincident with the results 
shown in Fig. [|. An important finding found here is that the boundary between 
the regions of effective and non-effective TPCF is located around Srj m 2.0. For 
relative momenta the boundary is located around AQ w 50 McV/c. Defining 
the region of effective TPCF in Sr] enabled us to directly study their intermit- 
tent behavior. Using the strip integral correlations we have driven a relation 
for the TPCF dependence on the phase-space pseudo-rapidity £77 in regard to 
the intermittent behavior. The dependence of logC^ m (A(5) on log AQ for the 
region 1 < AQ < 50 MeV/c supports our conclusion that TPCF show an inter- 
mittent behavior and probably add considerable values to the conventional FM. 
The dependence of the intermittence slopes <j) q on the orders q can confirm the 
conclusion concerning the intermittent behavior. To quantitatively estimate the 
contributions of TPCF we first studied the dependence of g-order C^™ on the 
second-order one. Here we noticed that the characterized behavior of intermit- 
tency, namely the ratios of intermittence slopes 4> q /(f>2 increase with increasing 
the orders q. These ratios are compared with results given in pi] . It is 
found that on the one hand the values supposed to be added by TPCF increase 
with increasing the orders. On the other hand the differences between their 
successive values decrease with increasing orders. 
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